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Abstract 

We rigorously establish some exact properties of reflection symmet- 
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ric spin systems with antiferromagnetic crossing bonds: At least one 

ground state has total spin zero and a positive semidefinite coefficient 

matrix. The crossing bonds obey an ice rule. This augments some 

previous results which were limited to bipartite spin systems and is of 

particular interest for frustrated spin systems. 
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1 Introduction 

Total spin is often a useful quantum number to classify energy eigenstates of 
spin systems. An example is the antiferromagnetic Heisenberg Hamiltonian on 
a bipartite lattice, whose energy levels plotted versus total spin form towers of 
states. The spin-zero tower extends furthest down the energy scale, the spin- 
one tower has the next higher base, and so on, all the way up the spin ladder: 
E(S + 1) > E(S), where E(S) denotes the lowest energy eigenvalue for total 
spin S ||. The ground state, in particular, has total spin zero; it is a singlet. 
This fact had been suspected for a long time, but the first rigorous proof was 
probably given by Marshall |l] for a one-dimensional antiferromagnetic chain 
with an even number of sites, each with intrinsic spin- 1/2 and with periodic 
boundary conditions. This system is bipartite, it can be split into two subsys- 
tems, each of which contains only every other site, so that all antiferromagnet 
bonds are between these subsystems. Marshall bases his proof on a theorem 
that he attributes to Peierls: Any ground state of the system, expanded in 
terms of S ,( - 3 - ) -eigenstates has coefficients with alternating signs that depend 
on the S^ -eigenvalue of one of the subsystems. After a canonical transforma- 
tion, consisting of a rotation of one of the subsystems by it around the 2-axis 
in spin space, the theorem simply states that all coefficients of a ground state 
can be chosen to be positive. To show that this implies zero total spin, Mar- 
shall works in a subspace with S^ -eigenvalue M = and uses translation 
invariance. His argument easily generalizes to higher dimensions and higher 
intrinsic spin. Lieb, Schultz and Mattis |J point out that translational invari- 
ance is not really necessary, only reflection symmetry is needed to relate the 
two subsystems, and the ground state is unique in the connected case. Lieb 
and Mattis ultimately remove the requirement of translation invariance 
or reflection symmetry and apply the M-subspace method to classify excited 



states. Like Peierls they use a Perron-Frobenius type argument to prove that 
in the S^-basis the ground state wave function for the connected case is a 
positive vector and it is unique. Comparing this wave function with the posi- 
tive wave function of a simple soluble model in an appropriate M-subspace || 
they conclude that the ground state has total spin S = \Sa — Sb\, where 
Sa and Sb are the maximum possible spins of the two subsystems. (In the 
antiferromagnetic case Sa = Sb and the ground state has total spin zero.) In 
the present article we reintroduce reflection symmetry, but for other reasons: 
we want to exploit methods and ideas of "reflection positivity" (see || and 
references therein.) We do not require bipartiteness. The main application is 
to frustrated spin systems similar to the pyrochlore lattices discussed in f7|. 

2 Reflection symmetric spin system 

We would like to consider a spin system that consists of two subsystems that 
are mirror images of one another, except for a rotation by n around the 2-axis 
in spin-space, and that has antiferromagnetic crossing bonds between corre- 
sponding sets of sites of the two subsystems. The spin Hamiltonian is 

H = H L + H R + H C , (1) 

and it acts on a tensor product of two identical copies of a Hilbert space that 
carries a representation of SU(2). "Hl = Hr" in the sense that Hl = h ® 1 
and Hr — 1 (g) h 1 where the tilde shall henceforth denote the rotation by n 
around the 2-axis in spin-space. We make no further assumptions about the 
nature of Hl and Hr, in particular we do not assume that these subsystems 
are antiferromagnetic. The crossing bonds are of anti-ferromagnetic type in 
the sense that H c = Y,a$a- Sa', with S A = E; 6 a ji St and S A > = T,i'eA'Ji Si', 
where A is a set of sites in the left subsystem, A' is the corresponding set 



of sites in the right subsystem, and ji are real coefficients. The intrinsic 
spins Si are arbitrary and can vary from site to site, as long as the whole 
system is reflection symmetric. We shall state explicitly when we make further 
assumptions, e.g., that the whole system is invariant under spin-rotations. 




Figure 1: Some possible crossing bonds. 



Any state of the system can be expanded in terms of a square matrix c, 



■0 = E C «/^a <g> ^/3> 
a,0 



(2) 



where {ip a } is a basis of S^-eigenstates. (The indices a, (3 may contain 
additional non-spin quantum numbers, as needed, and the tilde on the second 
tensor factor denotes the spin rotation.) We shall assume that the state is 
normalized: (ip\ip) = tr cc^ = 1. The energy expectation in terms of c is a 
matrix expression 

_ 3 
= tr cc ] h + tr(c t c) T /i - E E tr ^S^c{Sf)\ (3) 
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here (h) a/3 = {i}o\H L \i>p) = (^ a \H R \^p), (S^)^ = (ip a \ Eie^^sf %)> and 
we have used (S^ ) T = —(S^y. (For a = 1,3 the minus sign comes from 
the spin rotation, for a = 2 it comes from complex conjugation. This can be 
seen by writing S^ and S^ in terms of the real matrices S + and S~ .) Note, 
that we do not assume (h) a p to be real or symmetric, otherwise the following 
considerations would simplify considerably 0. 



We see, by inspection, that the energy expectation value remains un- 
changed if we replace c by its transpose c T , and, by linearity, if we replace it 
by c + c T or c — c T . So, if c corresponds to a ground state, then we might 
as well assume for convenience that c is either symmetric or antisymmetric. 
Note, that in either case we have (cr) t = cl, where cl = vcrf and cr = vdc. 
(Proof: (c 2 R ) T = (c*c) T = cc 1 = c 2 L , if c T = ±c; now take the unique square 
root of this.) Using this we see that the first two terms in the energy expecta- 
tion equal 2 tr c 2 L h and thus depend on c only through the positive semidefinite 
matrix cl- With the help of a trace inequality we will show that the third 
term does not increase if we replace c by the positive semidefinite matrix cl- 

3 Trace inequality 

For any square matrices c, M, N it is true that || 

| tr c ] McN ] | < - (tr c L Mc L M ] + tr c R Nc R N^ , (4) 

where c R = \fdc, cl = vcrf are the unique square roots of the positive 
semidefinite matrices (?c and cc^ . For the convenience of the reader we shall 
repeat the proof here: By the polar decomposition theorem c = ucr with a 
unitary matrix u and (ucru^) 2 = urfcv) = cc^ = c 2 L , so by the uniqueness 
of the square root ucrU' = Cl- Similarly, for any function / on the non- 
negative real line uf(c R )u^ = /(cc), and in particular u-^/cr = J~clU and thus 
c = ^/cEuy/c^. Let P = vly/cEMy/cEu and Q = ^/crN ] \Jcr~, then 

\tic ] McN ] \ = \ti PQ\ < -(trPPt + trQQt) (5) 



- (tr c L Mc L M^ + tr c R Nc R N^ , 



where the inequality is simply the geometric arithmetic mean inequality for 
matrices 

\ti PQ\ HE^'Qiil < ^E(l^l 2 + \Qji\ 2 ) = i(trPP+ + trgQt). 

4 Existence of a positive ground state 

Consider any ground state of the system with coefficient matrix c = ±c T and 
apply the trace inequality to the terms in {%fj\Hc\ ip): 

-trc^c^V > -i (tr c L S^c L (S^ + trc R S^c R (S^) , 

but (c R S^c R (S^) T = ((S^) T yc L (S^) T c L = (S^c L (S^)c L , so in fact 

-tr C t^c(4 a V > -tr c L Sfc L {S { ^. 

Since the normalization of the state and the other terms in (0) are unchanged 
if we replace c by cl = vcrf, and because we have assumed that c is the 
coefficient matrix of a ground state, it follows that the positive semidefinite 
matrix cl must also be the coefficient matrix of a ground state. 

5 Overlap with canonical spin zero state 

Consider the (not normalized) canonical state with coefficient matrix given 
by the identity matrix in a basis of S^^-eigenstates of either subsystem 

j 

- = 1212 12 ^(j,m,k) ® $(j,™,k>) (6) 

k,k' j m=~j 

j 

= 1212 12 i>{3,m,k)®{-y~ m %-m,k>y 

k,k' j rh=—j 



The states are labeled by the usual spin quantum numbers j, m and an addi- 
tional symbolic quantum number k to lift remaining ambiguities. The state 
S has total spin zero because of the spin rotation in the right subsystem: 
Its St ot -eigenvalue is zero and acting with either S^ t or S^ ot on it gives zero. 
The overlap of any state with coefficient matrix c with the canonical state 
S is simply the trace of c. In the previous section we found that the re- 
flection symmetric spin system necessarily has a ground state with positive 
semidefmite, non-zero coefficient matrix, which, by definition, has a (non-zero) 
positive trace. Since the trace is proportional to the overlap with the canon- 
ical spin-zero state, we have now shown that there is always a ground state 
that contains a spin-zero part. Provided that total spin is a good quantum 
number, we can conclude further that our system always has a ground state 
with total spin zero, i.e., a singlet. 

6 Projection onto spin zero 

Consider any state ip — Y^Captfa <8> i>/3 with positive semidefmite c = \c\. We 
have seen that this implies that ip has a spin-zero component. If total spin is 
a good quantum number it is interesting to ask what happens to c when we 
project ip onto its spin zero part 

V° = £c> a ®^. (7) 

We shall show that the coefficient matrix c° of ip° is a partial trace of c and thus 
still positive semidefmite. A convenient parametrisation of the S 3 eigenstates 
ip a for this task is, as before, a = (j, m, k), where k labels spin-j multiplets [j] k 
in the decomposition of the Hilbert space of one subsystem into components 
of total spin. Note that [j] k <g> \j'] k , = [j +/]©... [\j - j% so [j] k <g> \j'\ k , 
contains a spin zero subspace only if j = j', and for each k, k' that subspace 
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is unique and generated by the normalized spin zero state 

_i j 
6fe,fc' = (27 + l) 2 Yl Tpti^k) ® ipv,m,k')- (8) 

rh=-j 

(Recall that ip(j ; m,k') — (~y~ m i ) (j-m,k') is the rotation of ip(j,m,k') by n around 
the 2-axis in spin space.) The projection of ip onto spin zero is thus amounts 
to replacing c with c°, where 

if j '• ^ j' or m ^ m' 
j 

~ \ / y ^'(j,'in,k)(j,fa,k') 61S6. 



C (j,m,k)(j',m',k') — \ N ^ , (9) 

2j + l .^.' 

J m=—j 



(N is a overall normalisation constant, independent of j, m, k.) Let us now 
show that this partial trace preserves positivity, i.e., (v , c° v) > for any 
vector t> = (f(j )m) fc)) of complex numbers. If we decompose v into a sum of 
vectors Vj m with definite j, m and use (|J|), we see 

(v,c°v) = £K m ,c% jm ) = £ KLc<J > 0, (10) 

j',m j,m,m 

where the cu?^ are new vectors with components uTJfhk) = v (j,m,k), indepen- 
dent of rh. Every term in the last sum is non-negative because c is positive 
semidefinite by assumption. This result implies in particular that a reflection 
symmetric spin system always has a ground state with total spin zero and 
positive semidefinite coefficient matrix - provided that total spin is a good 
quantum number. 

7 Ice rule for crossing bonds 

The expectation of the third spin component of the sites involved in each 
crossing bond B, weighted by their coefficients jj, vanishes for any ground 

state ipo, 

(HEjM 3) + s¥ } Mo) = o, (ii) 

ieB 



provided that either the left and right subsystems are invariant under the 
spin rotation, h = h, or that their matrix elements are real (the latter is 
equivalent to the assumption h = h T , since we know that h = W or otherwise 



the whole spin Hamiltonian would not be Hermitean). By symmetry (II 
will also be true for the first spin component and, if we are dealing with a 
spin Hamiltonian that is invariant under spin rotations, it is also true for the 
second spin component. For ground states with symmetric or antisymmetric 
coefficient matrix we automatically have (s\ + s\/ ) = for any pair of sites 
% and i', so in that case (|TTD is trivial. 

For the proof we introduce a real parameter b in the spin Hamiltonian: 
H (b) = H — b(S B + S B ,) + b 2 /2, where B is one of the sets of sites involved 
in the crossing bonds of the original Hamiltonian H. Let E b be the ground 
state energy of H(b) and E the ground state energy of H. One can show that 
Eb > Eq and (|TTD follows then by a variational argument: 



(i> \H(b)\ip ) >E b >E = y> \H\iJ>o), (12) 

or, (ipolbisP + S { B ))\ip ) < b 2 /2, which implies (0). Note, that we did not 
make any assumptions about the symmetry or antisymmetry of the coefficient 
matrix of ipo here. 

Sketch of the proof of E b > E (see also |, 0): H(b) = H L (b) + H R {b) + 
H c (b) + b 2 /A with H L>R (b) = H LyR - b/2 ■ S^ and H c {b) equal to H c except 
for the term S B ■ S B , , which is replaced by (S B — b/2) ■ (S B , — b/2). If we now 
write the ground state energy expectation of H (b) as a matrix expression like 
(0) and apply the trace inequality to it, we will find an equal or lower energy 
expectation not of H(b), but rather of H: The trace inequality effectively 
removes the parameter b from the Hamiltonian. By the variational principle 
the true ground state energy of H is even lower and we conclude that E b > 
E . Role of the technical assumptions mentioned above: If h = h T , then 



the transpose in the second term in (|3|) vanishes, the matrix expression is 
symmetric in cl and cr (except for the sign of the parameter b), and the trace 
inequality gives (H(b)) c > \ {{H) CL + {H) CR \. If h = h, then we should drop 
the spin rotation on the second term of the analog of expression (0) for ?/v 
The matrix expression for (H(b)) is then symmetric in c and c T and we may 
assume c = ±c T to prove Ej, > Eq. The calculation is similar to the one in 
section |4j. Note, that c = ±c T only enters the proof of Eb > Eq, we still do 
not need to assume that the coefficient matrix of i/jq in (|ll|) has that property. 
The preferred configurations of four spins with antiferromagnetic crossed 
bonds in a classical Ising system are very similar to the configurations of the 
four hydrogen atoms that surround each oxygen atom in ice: There are always 
two hydrogen atoms close and two further away from each oxygen atom, and 
there are always two spins "up" and two "down", i.e. M = 0, in the Ising 
system. Equation (|TT|) is a (generalized) quantum mechanical version of this - 
that is why we use the term "ice rule" . This phrase is also used in the context 
of ferromagnetic pyrochlore with Ising anisotropy ( "spin ice" ) f| and we hope 
that does not cause confusion. 

8 Discussion 

We would like to discuss similarities between our method and previous work, 
in particular the approach of || for the bipartite antiferro magnet: There, 
the spin Hamiltonian splits into two parts H = Hq + H\. The expectation 
value of H Q with respect to a state ip = J2 fa<fia, expanded in an appropriate 
basis {(p a }, depends only on \f a \, and the expectation of Hi does not increase 
under the transformation f a -^ \f a \. The variational principle then implies 
that there must be a ground state with only non- negative coefficients \f a \. The 
present setup is very similar, except that we use coefficient matrices (c a p) to 
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expand states, since we work on a tensor product of Hilbert spaces. In our 
case the expectation value of H = Hl + Hr depends only on c via the positive 
matrices cl and cr, and the expectation value of He increases if we "replace" 
c by these positive matrices. The similarity is even more apparent if h has real 
matrix elements: In that case we may assume that c is diagonalisable and its 
eigenvalues play the role of the coefficients f a . The spin of a positive ground 
state is established in all cases from the overlap with a state of known spin that 
is also positive. In a system with sufficient symmetry we can, however, also use 
the "ice rule" to prove that all ground states have total spin zero j7|. (E.g., in 
a system with constant coefficients ji and enough translational invariance, so 
that every spin can be considered to be involved in a crossing bond and thus 

(3) 

in an ice rule, we would conclude that all ground states have S^l = and, 
assuming rotational invariance in spin space, S tot = 0.) It is not clear, if M- 
subspace methods can be used in the present setting to get information about 
excited states. An important point in the our work is that we consider not 
only antiferromagnetic bonds between single sites but also bonds between sets 
of sites. This frees us from the requirement of bipartiteness and even allows 
some ferromagnetic crossing bonds, for example in (si — S2)(si' — s 2 '). There 
is no doubt that the scheme can be further generalized, e.g., to other groups 
or more abstract "crossing bonds" . In the present form the most interesting 
applications are in the field of frustrated spin systems M . 

We did not address the question of the degeneracy of ground states. Clas- 
sically a characteristic feature of frustrated systems is their large ground state 
degeneracy. For frustrated quantum spin systems this is an important open 
problem. 

We would like to thank Roderich Moessner for inspiring discussions. 
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